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by 
JOHN DENNIS BOWEN 


Submitted to the Department of Ocean Engineering in partial 
fulfillment of the requirements for the degrees of Naval 
Engineer and Master of Science in Naval Architecture and 

Marine Engineering. 


ABSTRACT 


The specific impetus for this work was a conceptual 
design of a submarine using the toroid as the pressure hull. 
The designers could not find a ready body of knowledge to 
obtain scantlings for thier pressure hull. 

This work began with a review of efforts to solve 
complete toroidal] structures. Several works were found which 
addressed general shells and extended into partial toroids, 
but. the solution of a complete toroid was not found to be a 
common exercise. Some of the these works are briefly 
reviewed. An attempt was then made to solve for the 
displacements in a thin walled circular toroid using the 
energy method. Several problems were identified associated 
with the structure geometry which make the solution for the 
complete toroid difficult. In addition. the functional used 
for the energy method needs to be more complex than the 
Simple trigonometric or power series functionals used in this 
work. These two areas, geometry and functionals, are fertile 
areas for further study 
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CHAPTER 1 
INTRODUCTION 

The purpose of this work was to determine the response 
of a thin toroidal shell structure. The initial goal seemed 
achievable until the literature was reviewed. A simple and 
understandable solution was not available. The disparity in 
works was apparent in such basic elements as coordinate 
systems. The solutions were generally directed towards 
displacements or forces. The solution paths varied from 
complex algebra to finite difference method. The solutions 
were as varried as the number of individuals attempting them. 
lt. was impossible to expect any simple correlation among the 
few people who attempted to address this topic. 

The result of this work is the identification of two 
areas where more time and effort is required. 

A. The impact of geometric curvature on the solution. 

B. The impact of the shell geometry on the assumed 

functional used in the energy method of the solution. 


To be very specific, these two areas have precluded this 





effort from achieving the initial purpose of this year long 
cGDROrt. 

A word needs to be addressed towards assumptions that 
are made. The assumptions are fairly standard for a shell 
being analysed in the linear elastic range. The assumptions 
are: 

1. The shell is made of isotropic and homogenious 
material which obeys Hook’s Law. 

2. The thickness of the shell is constant. 

3. The thickness of the shell is small compared to 
the radii of curvature. 

4. A straight line normal to the middle surface 
before deformation remains straight and normal to the 
middle surface after deformation and retains its 
original length. 

In order to establish a consistent presentation, 
geometry will be discussed here and, unless otherwise noted, 
all representations of geometry will follow this coordinate 
system and the definitions. Solutions by other individuals 


will be translated into this reference system. 





DEFINITIONS 


a 


Fatro of E to r; a= > 


Since R > r to form a complete toroid, a will always 
be greater than 1. 
Circumference of the parallel circle: C=2n pe 
Radius of curvature in the © direction of the 
unloaded shell envelope, or radius of the circle 
which is rotated about the Z axis (at distance R) to 
form the toroid. | 
Radius to any point on the unloaded toroid measured 
perpendicularly to the Z axis. 

r" = R +r sinO =r Ca + sin). 
Radius of rotation about the Z axis of the circle of 
radius r used to form the toroid. 
Thickness of the plate 
Angle of rotation about the Z axis measured counter 
clockwise from the positive X axis in the X-Y plane. 
Anzle of rotation measured clockwise from the local 
perpendicular to the X-Y plane in the positive 2 
directon at distance R from the origin in any 
direction 6. See Figure 1. It sounds like a 


tote. Dut it is not. 





Figure 1 


Layout of Variables 





ede 





p Radius of curvature in the og direction. 


ANIE 


= + r 
sine 





Coordinate System: The standard right-handed XYZ coordinate 
System 1s the global reference system. Other 
coordinates Сес. Ө апа @) build on this basic 
reference system. 

Geometric Curvature: Curvature of the shell due solely to the 
geometry or curvature of the unloaded toroid. 

Load Curvature: Curvature of the shell caused by 
displacements in response to the applied load. This 
could also be looked upon as the change in curvature 
from the reference Cunloaded) curved surface. 

Axis of Symmetry: Axis about which the small circle is 
rotated to form the toroid - the Z axis. 

Plane of Symmetry: Plane which bisects the toroid - the X-Y 
plane 

Meridian: Intersection of any plane containing the Z axis 
with the toroid. On either side of the Z axis, the 
intersection results in a circle (for the unloaded 


toroid) with O as the variable. (See Figure 2) 





Parallel Circle: Intersection of any plane parallel to the 


X-Y plane with the toroid. (See Figure 2) 
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Figure 2 


Parallel Circle and Meridian 


Y 
| Parallel Cire le 
(5 f 


ин. 2 
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CHAPTER 2 
HISTORICAL VIEW AND THE PROBLEM 
Historical View: 

Upon starting to research the analysis of circular 
toroidal shells, the field was found to be not very large. 
Books on shell theory and analysis would neglect or only 
brietly mention the complete toroid. This could imply that 
the toroid was: 

C150 a simple extension of general shell theory or 

(2) much more complex than general shell theory. 

The latter seems to be the case. 

Senjanovic [1] gives a brief history of shell theory 
in his book. His Russian pride shows in some instances 
c"...shortcomings were eleminated by representatives of the 
Russian School..."), but the relative youthfulness of the 
field is obvious. Initial work in what is now called shell 
theory was by Н. Aron (1874) and A. Love (1888) a mere 
century ago. Sen janovic does cite Flugge ~ another work 


referenced here - as one of the developers of the classical 





Cor western) theory of shells. Yet even Flugge devotes only four 
pages [2] to the toroid. The following works are illustrative of 
the background on toroidal shells. 

W. Flugge: Flugge starts with a generalized shell and 

establishes force equlibrium. From this equlibrium 

he solves differential equations to obtain solutions 

for forces and displacements. His discussion of the 

toroid points out some of the problems inherent in 

the full toroidal solution. 

L. Sobel: Sobel [3] was an understudy of Flugge. 

His doctorial dissertation starts with the same 

generalized shell theory but is slanted towards the 

buckling analysis of the toroid. His solution was 

worked on a computer and is difficult to follow. 

I. Senjanovic: Senjanovic starts with a force 

ballance on a differential element in a localy 

orthogonal coordinate system, then uses several 

methods to solve them. Most notably, he employed 

complex variables and some computer solutions circa 

1970. 


©. Timoshenko: Timoshenko [4] starts with a plate 
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and develops into a generalized shell. His attention 
to the toroid is of the same order as Flugge. This 
writer found Timoshenko's general development of 
theory easier to follow than others and hence the 
references to Timoshenko in this area may be more 
frequent than to other authors listed here. 

E. Tsui: Tsui [5] starts with the constitutive 
relationships and develops them into forces. However 
for the toridal shape his solutions are on globally 
orientated forces and displacements. Tsui intended 
his work more as a handbook than a textbook. His 
solutions are tables of values from computer analysis 
of toroidal shells. He does provide however, a 
consistant development, though shorter and in less 
detail than some others, for various shells of 


revolution. 





The Problem: 

The specific problem to be addressed in this thesis is 
the response of a toroidal shell to hydrostatic presure. 
More specifics on the loading will be addressed in chapter 
three. This is to set the stage for the problem itself 
and lay out the difficult areas. 

The goal of any structural solution is to take into 
account the loading, geometry of the structure, and the 
material in the structure and determine the response of the 
entire structure. The response is then compared to a desired 
standard or a set of failure criteria and adequacy of the 
design is determined. This approach is generally slanted 
towards displacements of the structure since only through 
displacements are the applied loads distributed and 
equilibrated. 

In any structure, the various elements are in 
communication with one another. The beam under a simple 
tensile axial load is the easiest example, but a plate or 
shell must also satisfy this connectivity. Unfortunately, 
the complexity of this connectivity or communication among 


elements increases rapidly. In the axially loaded (tensile?) 





beam, only the axial displacement «say the X direction) is 
important. Load the same beam laterally and now the X and Z 
directions are required to describe the beams response. 
Going to a plate now requires the X, Y, and Z directions. 
The shell now takes the plate out of the convenient XY plane 
for a reference. 

With appropriate coordinate system selection, some very 
practical shells can be handled easily. The cylinder and the 
sphere are classical examples. The transition to a toroid 
would - by simple implication - be very easy. Start with a 
cylinder of radius r and length of 2nR and bend it around 
upon itself. Discounting a few wrinkles on the inner 
sections of the cylinder, we now have a toroid. The wrinkles 
however carry a much deeper implication. 

The structure’s ability to distribute the applied load 
changes drastically. The complete toroid has no boundary 
conditions comparable to the cylinder’s boundary conditions 
at X=0 and X=L. These are replaced by connectivity Cor 
continuity»? requirements. The analysis of toroidal segments 
is facilited by the ability to insert boundary conditions. 


But even this facility can be of limited use if too much of 





the structure is included in the analysis. The simple 
thin walled analysis on a cylinder for a force balance 
becomes an algebraic exercise of no small feat for the 
toroid. This will be shown in chapter five.) 

It should be noted that all of the individuals above 
st.ated or implied some difficulty in obtaining a solution for 
a complete toroid. The most blatant was Tsui [6] when he 
states that. in problems involving the complete toroid, the 
crowns € Ө = O and 180° - see FIG 1> can not be solved 
due to singularities. Flugge [7] gives a solution for the 
forces but states that they cannot be used "... in the 
vicinity of the top and bottom circles" CTsui’s crowns) 
"without additional bending...". He then provides a 
displacement expression containing terms which blow up at 
Өө 0° and 180°. 

With all of this going against an easy solution, two 
courses of action were undertaken: 

A. The singularities at the crowns could be avoided 
by using some expression which avoids the offensive 


term С > (see Chapter 60. This was avoided by 





sing 


referencing the differential element to: 


= 17 - 





г = Е + гап е: г са + ѕіп Ө > 
The biblical addage still holds: If thine eye offends 
thee, pluck it out. 
B. The method of solution selected would be one to 
avoid further problems as pointed out by Flugge [8]. 
The energy method would avoid messy mathematical 
stumbling blocks in the solution, since the form of 
the solution is assumed going into the energy method. 
This was almost true as shall be pointed out latter. 
Additionally, this writer has not seen this attempted 
by others. Therefore, when the line forms to Say 
that. this method will not work, it should be a very 


short line indeed. 
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CHAPTER 3 


CONSTITUTIVE RELATIONSHIPS 


An expression. or series of expressions, relating the 
material motion to properties such as stress and strain Cand 
further to forces and moments) is required. This follows the 
definition of constitutive equation of the material as given 
by Fung [9]. "...mathematical expression of the mechanical 
property of a material.” 

First we must define our displacements Csee Figure 3). 
Three displacements will be used: 

u :Displacement in the positive © direction 

v :Displacement in the positive € direction 

w :Displacement in the direction of the positive 
Cie. outward pointing?) normal at the point of 
interest. 

Rotations will be defined as follows Csee Figure 3): 

dg : Rotation of a unit normal about the tangent to 
the meridian in the positive u direction: that is. 


rotation about the positive © axis. 


EN. ue 





Figure 3 


Dispalcements and Rotations 


Displacements 





Retations 
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On ` Rotation of a unit normal about the tangent to 
the parallel circle in the positive v direction; 
that is, about the local positive q axis. 

o. : Kotation about the positive w direction. 

The loading will impact greatly on the final outcome. 
For this problem. the loading will be assumed to be a 
constant hydrostatic pressure. This has several 
implications: 

A. Hydrostatic loading always presents a loading 
normal to the surface. 

B. Because of the global nature of the hydrostatic 
loading, anv benifit to be garnered from symmetric 
conditions would be applicable to this problem. 

The second aspect of the loading, axisymmetric loading, 
results in the number of unknowns being reduced dramatically. 
Both Timoshenko [10] and Flugge [11] state that the total 
deformation can be represented by only two displacements; one 
normal to the surface and one along the meridian in the 
positive © direction. Flugge also states [12] that that 


stresses are independant of @ and in fact any derivatives 


with respect to @ are zero. Timoshenko 113] agrees. but. 
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en zderivative going to zero is implied, not stated as 
bluntly as in Flugge. 
To take this one step further, a lot of things will go 


to zero because of axisymmetric loading. Sheer stress and 


strain are cross coordinate terms, ie. 


Be 2 
Dela, = $5 € рж] ) 
EN SG SO 
| 2577 _ E 
But any Só 0. SO LE] og = L*] 66 = Q. 


Because of the type of loading described, v and © are not 
required for the description of motion. v is not required 
because of the independance of the motion on @, and O, is not 
required because it depends on a partial of 4. Allof this 
is the result of axisymmetric loading. 
STRAINS 
Midplane strains are as follows: 

Eg : ls the strain in the positive 6 direction. £o 

is made up of two parts (see Figure 4): 

1. £g due to u 

At two points on the meridian; P and Q 
ou 


usu u, u + 35 de 


SA + BU 16 — u) = 
1 га 50 


D £o due to w 


x 
ФЕ 
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eo 


Figure 4 


Midplane Strains 
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Length of a segment is equal to the. radius 
times the angular rotation. The original 
length is r dO. At point P the radius 
increases by w to Сг + w) so the new length 
would be (г + и) 90. Similarly at point Q the 
radius increases but is now 

(TW o де). The new length would be 

Cr +w + 2146096. Averaging the two lengths 
at P and Q gives an average Al = wdo + > Sido. 


Dividing by the original l = r dO gives: 


с = - ім + 2 ŚW do] Neglecting second 
2 60 
order terms yields: = = 1 


Total eg is therefore: 
со = т 152 + м) 

Ey: Is the strain in the positive @ direction 

Csee Figure 4). Since the body is 

symmetrically loaded, all 6 s have the same 

displacements. Therefore Ep can be calculated 

by calculating the change in the circumference 


of the parallel circle at 6 = constant. Since 


the circumference is related to the radius, 
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the radius a will be used to determine 
straln Eg: 

C=2 nr“ 
At point P, displacements w and u contribute 
to the Ar” as follows: 


Ж | 
Аг = u cos@ + w sin® 


AC = 2 т Аг” 


STAG M : 
Es FA ms lu cosO + м sinO] 
ROTATIONS 
Op is determined from two points, P and Q, on a meridian 


seperated by displacement u in the positive € direction (see 
Figure 5). The unit normal vector at P must rotate 
through do to go from P to Q due to displacement u, ie. 


г ДӨ = u 


If w was constant from P to Q, no additional rotations would 
be required since the surface displaces uniformly through dé. 
If = z 0. then a rotation of the local normal at 0 will be 


required to maintain its normal relationship. This rotation 


will be: 


N 


о, = — ŚMd6/(r d6)= - IM 
68 róo 
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Figure 5 


Rotation Og 


w w- Wo 
wy € = / 
Q 


Jo 
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This makes the total Oo 


бм 


CES 1 (u = 
er 68 
Because of the symmetric loading and no subsequent 


variations in the displacements with 6, no v displacements or 


ów 


exists. However a ig Can exist (see Figure 6). 
[9 
eo 


Due to symmetry in loading, Iwal is constant in @. 
However. for o, to change direction. an additional rotation 


Ф 


perpendicular to Op must exist. This is the rotation Og. 


сөз - Op cosé do 
Figure Ó gives a graphical presentation of Qo. 
Timoshenko [14] gives a similar result but achieves this as 
a Ec result of obtaining curvature. 
BORDZCURVATURES 
Load curvatures are as follows: 

ха: Change in curvature along the positive © axis 

due to the applied load. 

Ха: Change in curvature along the positive Фф axis 


due to the applied load. 





Figure 6 


Rotation Gg 





4 (0 = de, sın(d6) = аж, de 





/ 


+ 
4 
+ 





-4W lose 
— (wg 4%) lose 
= de Cose de 
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It 1s significant, to this writer at least, to seperate 
these curvatures from the geometric curvature defined in 
Chapter 1. The geometric curvature gives one an indication 
of how difficult the problem will be [15] and defines the 
initial or unloaded reference surface for further work. The 
load curvature is the sole result of the structure’s response 
to the load applied and directly related to the displacements 
use to equilibrate the load. 

Load curvature along the © direction is the change in 
rotation about the 6 axis b divided by the original 
length (rd6). 

"o = w S Op 46 


0005 


Recalling: og = 2 (u - ŚW) 
Г 


1 Фи 52м. 
СТИ Са - =) 
Ө 127506 66? 


Timoshenko [16] gives a description of this and x is 
Similarly: 


Og дехсг“ de) 


= 
! 


cos 6 ÓW 
ONZE i (CU = ) 
r2(a + sin0) oU 


Ve can now describe all the pertinent engineering 
parameters in terms of displacements. Adding to this list 


the general expression for planer stress, we get the 





following: 


де 


9% 


This 


E 
pp 
E 
1” 
Е 
lear 





2 





2 





2 


almost 


Stress, strain, 


expressed in terms of 


Step is to proceed to 


follows: 


C£g + DE 4? 


E ES + w} + ——____ [u cos® + м sin®)) 
rca + sin® 

qu [u cosO * w sinO) + ES (eu + wJ> 

rca + sin® r 


completes our constitutive relationships. 


rotation, 


h 

N =,“ 
=h 

2 

h 
"EU 
V -h 
2 


and local curvature can be 


the displacements u and w. The last 


normal forces and moments. This is as 


o. dz ; i = Ө, e 
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CHAPTER 4 


EQULIBRIUM OF FORCES 


The next aspect of the solution to be presented will be 
that of establishing equlibrium conditions for the 
differential element. The differential element is shown in 
Figure 7. Geometry rears it’s ugly head in that the 
element. needs two reference points to describe it’s surface. 
The first is the axis of rotation of the shell in the © 
direction. The second is the origin of the global coordinate 
system. The first axis can be referenced to the global 
coordinate system by R and ©. 

Solutions for partial toroidal shapes all take advantage 
of the Known boundary conditions on the edges. When a 
complete toroid is analyzed, the boundary conditions double 
back upon themselves and become internal vice external 
forces. In these instances, symmetric loading and 
deformations are invoked. Tsui [17] gives tabulated 


computer results for both solutions and the differences are 





Figure 7 


Differential Element 
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obvious. 
In the following discussion, and the associated figures, 
the following conventions will be used Csee Figure 8): 
Р, is a generalized force ie. a force or moment in 
the i direction G = Ө, 9, n>. The units for E are 
[Force/unit length). 
Р. іѕ “һе Ғогсе, іе. 
F = [Forceunit length)*[Llength] = [Force] 
o: Ny Normal forces in the (5 direction acting on 
the side of the element on which € is constant. 
Qo Sheer force acting on the ф = constant face of 
the element. Note: Due to symmetry of loading, 
Mo- Mo Moment to induce rotation in the с» 
direction. The moment is defined as positive with 
respect to the right hand rule as applied about the 
local positive <2) direction. 
Each force will be addressed graphically with three 
isometric diagrams and a geometric breakdown to the 


appropriate Е. The F for each force No: No . .ect.5 can 


then be summed for a force balance. 
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Force Figure 


No 8 
Ny 9 
Qe 10 
P Cload)> 11 
Mo 12 


МФ 13 





Figure 8 


Forces Defined 
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Figure 9 


No 


> боф 


Ве в = 


ме 


d¢ 





Contribution to po doe to Ng: 


? > ON > > 

F,: -Na Sinf + CN, +? — 940» cosa sinß 
Ah © Ө бө 

е ‚ 6Ng А : 

Fa: -Мл сос) + (N, + ——d@) cosa sinf3 
© © © 58 


> > No > 
Kae N. sina ~ (N, + —d® sina 
© ө е 
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Figure 10 


No 





Contribution to F doe to №: 


Еф: “Ny + № cosa cos 
Po: No sina cosß - № cos(3 sina 


p No sin - No sinf3 
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Contribution to Е, doe to Qo: 


Fo Qo cosa sin + Qo cosa sinß 
Po: On sina cos - 96 cosa sinß 


Fn Qo cosa cosß + Qo cosa cosh 








Figure 12 





Р 
я 
= | 
P e 
бо 
1 
a = Z ( 40 + o, + 5520) 
EGZ 4 - Cw o 
2 © S$ no 


NOTE: Z = P Area = P [CR ч г віпӨоа4Ф1 (г де) 


Contribution to Е, doe to 2: 


Fo: - Z cosa sin? 
Fo: — 2 cosß sina 


En: -  cos(3 cosa 
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Figure 13 





Contribution to F doe to Mo: 


> 


, , ФМ ф , > 
Ғы: Mo cos - Mo + Py de> cosf3 cosa 


> 


Р óMg 
Fu: -(M, + ——d® cos sina 


; , бм 
Е: - Mg sin@ + (Mg + O 
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Figure 14 
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Contribution to БЕ doe to М 


ф 
> > ФМ > > 
ЕНЕ M, GOS cosa — CH, + ur dO» cos(3 cosa 
eo Ф = 
; ФМ 
Fa: - М, сос) sina + CM, + — qe» cos(3 sina 
е Ф 97% 
, > ФМ › > 
FL: - M. cosa sing + (M, + —£ de> cosa sinß 
n w Ф Se 


lo 





Relationships between the N and N ¿must now be 


established. 
-— Ñ ж 
SN , ON, ж 
No + OEN өз ea oe do 
50 50 50 
; ; ж ÓN ÓN xk 
= Nardo + No з»аф + с- Фг”аф + <- 2 «ФГ ag 
ге 60 60 59 
> ON жк 
= “No Go S quer dd, neglecting terms wither as 
50 
being relatively small 
Similarly: 
No, = № Ido 


Qo = Q o T” dọ 


DO, + E rag 
Мо = М г dé 
SM M. 
REŻ s Mo —9» r" ae 
SO 
do 





CHAPTER 5 
ENERGY METHOD 
Any attempted solution must maintain some sort of 
equlibrium. In the energy method, the equlibrium comes from 
minimizing the energy from the external load and the internal 
energy of the shell. For the external energy, a surface 
integral of the hydrostatic load and the normal displacements 
gives this value. The internal energy is determined by 
integrating stress times strain throughout the shell volume. 
The external energy is subtracted from the internal energy 
and the resulting functional is generally refered to as Il, 
the total potential energy of the structure. With V 
representing the internal energy and Y representing the 
external energy, the functional looks like: 
Hav- W 
If the value of N is at a minimum value, the structure 
will be in equlibrium. This is the principle of total 
potential energy [18]. This work will not address buckling 


or post buckling behavior, only the linear elastic region. 





For this situation, equlibrium is achieved when: 


бу óW 
og Фа Фа 


V V v 


= 0 


where q, is any generalized parameter used to determine V and 
E. 

Unfortunately. the parameters used to determine V and W 
are u and w — which is what we are attempting to solve. To 
get around this difficulity, a solution will be assumed. 
Then an energy balance will be conducted. 

The assumption of a solution is not without risk. 

First. the solution must meet Known natural boundary 
conditions. Then the accuracy of the solution is limited by 
the closeness otf the assumed solution to the actual Cand 
unknown) solution. 

In the case of the toroid, there are no “boundary 
conditions” as was pointed out in chapter four. In place of 
the boundary conditions, we have conditions of compatability 
to limit the solution. As stated earlier <119] апа [20]1», 
symmetry of deformation can be assumed for the hydrostatic 


loading. The plane of symmetry is the XY plane. 
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Eq 


Eq 


Eq 


Eq 


Eq 


Eq 


be 


To determin 


Ve= 
= 
re 
У = > e 
qu No Eg = 
Similarly: 
2: М) 2 = 
зи Ме Хо + 
d: Mo X y, = K 
Recall from 
| уі 
9.06: te => r 
1 
TB a = — 
© r2 


The unknown 


made that ua 


e the internal energy we proceed as follows: 
& j = Ө, ф 
2 45 


Recall that, due to the type of loading, 


cross terms €. Æ j) equal zero. 





C + + + 
No £ Ny Ey Mo X p My x) ds 
| E h 
+ 5 = 
К Eg v Ey? Eo ; K ae: 
2 wa 2 
K (eg + 2 р Eo Es + Eg? 
Е 2 2 
K сє + 2 v Es Eo + ES 
E h? h? 
D Cza + x,» x ; Dang = K сә 
ө С СӨ 12<1-v2) 12 
K ch’, Cx + 2 van x, + 22d 
12 Ө ө “о Ф 
LES EIA xo 
12 ф p “O © 
chapter four: 
cou + w); Es Ера Пор Cu cos® + w sin® 
68 rca + sin® 
2 
оз - ŚW; ky = zy cu - 2% 
го SO r“Ca + sin@) бо 


s u and w still persist. The assumption will 


nd w are functions of Ө. A series solution 


will be assumed as follows: 
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мсб) = 2 w sinćmO>, limits of summation will be 
addressed later 
u(O = Z u & csincmer> 
" 60 
ш 2 


u m cosmo) 

This turns out to be not a very good assumption for the 
the series representation of the displacements. There are 
some inherent limitations; for instance sin€m®) is always 
zero at, @© = 0. While not in conflict with the conditions of 
compatability, this is an additional constraint on the 
solution displacements. If the actual displacements are not 
zero at O = 0, then there is a guaranteed error in the 
assumed solution. The question to be answered, as with all 
assumed solutions. is how good is the answer obtained 
compared to the answer that is required. More on this 
subject in Chapter 6. 


For ease of notation, let: 


sinCimo» 


Ф 


m 


? 


Ф = m coscmed 


m 


ó = -m* sincme> 
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Therefore: 


Eq 9: 1-2, 


! 
M 
T 

$ 


Eq 10: “Өө = зц 6. 
Boundary conditions for this problem are actually 
compatability conditions. At Ө « 90% and 270% the uce> 


displacement should be zero and Ô сисе»> should be zero. 
Ө 


ОП о сн фә = зм О сфу ау н д 


Since w, are merely coefficients now in the series solution 
they are unaffected by the derivative. To look at this in 
more detail: 

5 EUM costo) + We 2 созс205 + м. 3 созсз0» +. 


*w n cosCno» 


For conveince, let n = 10 for the trial solution. Аб В = 909 


and 270. cos(m6> = 0 for m = 1, 3, 5, . . and cos(mó> = + 1 
гет гс 2, 4. 6. . . . For the series to be equal to zero 
set w = О form = 2, 4, 6, . . . Similar reasoning gives 
Nut Ücporambw 2. 4. 6. .— CThis explains why the 


functional for u was selected as the derivative of the 
functional for w .> As shall be seen later, the method is 


smart enough to realize this limitation.) 


The elemental area CdS) must be investigated for the 
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toroid. dS = dl, dl dl, = r do =; dl, = (R + r sin® dọ 


ф 5 
= r dO (R + r <1һбӨ› аф 

As a check: Ss Jama СЕ + r sin0> de de 

= 4 n? Rr. This agrees with the CRC 


Standard Mathematical Tables [21] for the surface area of a 


toroid, therefore dS is correct. 





The ultimate objective is to get on = O. In this case 
4 
а = ч апа “м. So for each specific Uu that is u, — one 
gets an expressison such that: 
ОП - ФУ . óW no 
Su, Su, Su, 
Likewise a series of m expressions: 
oll - ФУ _ ów _ 0 
Ow ÓW, ÓW, 


For a nominal m = 10, this yields 20 equations. Fortunately, 
u and w are 20 unknowns. 

Now the internal energy can be formated. Equations 9 
and 10 are substituted into equations 5 through 8. These 
equations are then inserted into equations 1 through 4. This 
is combined with the expression for dS to give V, the 
internal energy. To facilitate the future derivatives which 
will be required, the subscripts т апа п will be used to 


indicate the seperate terms in a product, ie.: 


San 





ES Eg = = Cu +Em ọ > 


m m m 


(Zu $ +iw ¢) 


2 


All ot the above yields: 


„2% прай B 1 : 
IST SE: So [zu $9. + >м $ Cou pP + 2n 9) 


NE 1 , 
po eu ZWm?m Pla F SINOSCOSY >Ч Фф + sinó 29 P, > 


1 > 
Gam Sino cone 2u d$. + sino 2 wi? > 


X Ta F-sręgy(Cos© Xu $ + sinO У м ф > 


n^? 1 . +8 4% 1 a6 su 
+ q? -z C2u o. * >м $ > ADO + 2w $ 2 

ho (22 s > coso 4 $ 
+ (>) — (2uoó + ZW $ I)— (ud +2uo) 

12 г2 m m m m г2<а + sine пп пп 
созе , s 
2 DCU ŻW O DV 
r“Ca + sin®) а Sa 
с. 030 Su $ + Zw 6 >] Ca + sine) de 
r“Ca + sine фа А 1 


Next. carrying out the multiplications yields: 
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A н 1 = es oe oe oe 
жама n K го С 226u M. A р en рав T E Vie 
И "и e ae a 


+ 2v 22Cu u M cosO + uw фф. sinó 


r*Ca + sino) 

* Wu Ф 9, созе + “ww Ф 9. sino» 
-n cuit S dd u ó соз2е 
r^«a * sine»? ENS капы 


+ uw фф. созе ѕіпӨ + wu фф cos® sin® 


2 
+ Wow Ф 9, сіп “ӨО 


h* кор», .. . oe ee 
ma e, ека РЕ Б иса 


+ с 





ом, 999 


A %,% 


m N m m 


2 3 
+ С h T .2 v COSO sscy u 
12r Ca + sine Ше” 


ве » ve > 


т 





n? > cos^e 


> > oe > 
+ < Е złu u ę pb = U wę ОФ 
12% Ca + sinó», 





mn m m m mn m m 


- wu фф, * WW Ф Ф >] Ca + sin® do 


mn 


This is all well and good but actually oe is what we 
oq 


L 


really need. Again, for this problem а =ч and w and 


k k? 
k=1. 2, . . . 10. Applying term by term differentiation is 


straightforward but tricky and tedious. A couple of 


examples: 


- 50 - 





But. observe: 


= U d. d, * U, ФФ. 


ZN 
е 
3 
е 
= 
© 
qi 
© 
S x 
Мл 
1 


2 u, 6, 9, 

Since m and n can both serve as a counter index, the two 
terms may be combined. The key to this is that both 
functionals cp) have to be of the same order of the 


derivative. Observe: 


5 et > a 4. > ee > 
óu cu, u. e. $? a ци = e. E ur ФК e 
E 


+ На. СФ Ф, E Ф, > 
Doing all of the above and collection like terms gives the 


following two expressions: 


The partial of V with respect to u. 


OV Е 21 Е te 19 2 1? cose ee > > ee 
Su, ы, To о 
2 > > 
B 2 созге ф Ф, 


Ca * sino»? ™ 
h? "^" 2 v cose ? á 
+ PE Omk ta F Sinos fkm SINO + ФФ, с0502 


2 сов2Ө 


2» ee sine 
Са + sin® 


m 


27 р в. 
зам 2 (29,9, + 


2 cos® sin® ,’ 


+ © >) coso) + p, $ 
р Са + sino? К т 
h? : 1 2 1? cose ae > te > 
+ C 2<-2 фф - ------- (офф + ФФ)? 
2 cos*e 


= от $ „Ca + sine) de] 
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The partial of V with respect to м 


Фу : 27 ji 2 v = , 
— = nk u Г^'С2 ф $, * ———————$ Фф, sind + ф ф cos® 
ÓW, [Rum Lo gut us Ca * sino "X DE 
> 2 cos® sin® |,’ 
(a + sino? ™ * 
Е necs "2" 2 р coso Mh PS 
C77 Ce PP; Ca + Sindy PP E PP? 
2 cos*e E 
Б —RUe 200a ct sino» de 
Ca + sin® 2 


4 v sino 


+zw 72% (2 фф + 
ши. ^ "  €a + sine 


$ 1 


m 


2 sin*e 


+ 
Са + sino»? 


Ф, ф 


m 


het bur 2 Ł COSO sad Sm 
СС Cm RP? 


2 cos*e 


ŻĘ 
Ca + sino»? 


6.6 >>Ca + sin®) de] 


The external energy (W should be looked at next. The 
energy is the product of the load - hydrostatic pressure 
(P> in this case which is always normal to the surface - 
times the displacement along the line of action of the force 
- wCOD - in this case. 
We sf P w CO» dS 
= JO" gd" PZ w, 6, r СЕ *r sinO dO dé 


Substitute R = r a 


ОР? м. „= p Са + sine) d® 





And the partials are: 


ov = 2nP r? 52% % Са + sine de 
бм 

k 

бн Ern 

óu 


k 
The series expressions can now be rearranged: 
ón ФУ _ ów 
слана от 
ФУ _ ÓW ) 
а Sq 2 i к? "к 


l t 





This can be represented in matrix format as follows: 


Eq 11 ОУ = [DJ <u> + [E] «w s 0 = OF 
ou Su 
k k 
Where [D) represents the coefficients of u іп OV 
a ou 
k 
(see page 9). 
[E] represents the coefficients of W, in ey 
ou 
k 
(see page 9) 
al = 
{u) = u, 7 U,, U,, U,, uio 
ое = 
(Ww) = Wo 7 M.» MS, Ta, М/с 
Similarly in matrix format: 
Eq 12 V e [F) <u> + [G] <w} = £P) « OF 
ÓW ÓW 
k k 
Where [F] represents the coefficients of u іп ON 
» ÓW, 
(see page 10). 
ГО) represents the coefficients of w, in За 
“ч 


(see page 102 





{m one. OW а ONE. gow 


ÓW, ÓW, ÓW, Mio 


Solving equation 11 for {u} gives: 
{u} = -— [DI] (ЕЈ <w) 
Use this expression in equation 12: 
«IFJ EDI * [E] + [G]? <w = <P> 
Solve for tw): 
tw) = ([F] [DI] * [EJ + [G]> * <P> 
If the solution for W and u_ are correct, then these 
can be used to reconstruct wC®) and uc® and from these, 
using the constitutive relationships in chapter four, 
stresses in the © and 6 directions can be calculated. Roark 
[22] gives expected values which can be used for comparison. 
Some justification for Roark*s results should be 
provided since this is how the solution will be checked. He 
provides no specific reference for these answers in his book, 
but working backwards one can see that a simple thin walled 
approximation analysis was performed, similar to that of a 


cylinder. Roark provides: 


t a * sinO 


(See Figure 15) 


a = 





Figure 15 


Thin Walled Analysis 
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Од is a simple force balance in which the pressure times the 
internal cross sectional area equals the stress times the 
material area. 


PA = 0% A 


m 


Р n r?) = ФАО? 


7e 
Pr 
a с=з 
им 
A similar approach is used on de: The pressure acts on the 
area A, which in this case is (2 го<2 т КО. O6 acts at O and 


Ө + 180%. The thickness is constant at t but the 


circumference at € and 0 + 180? are different. 


CCG) = Ż nr” = 2 n r(la + sing 
The metal area at © is therefore ACO) = CCO) t. The force 
balance is: 
РА - о<Ө,2 С<Ө,2 1 жоСӨ,2 С<Ө,2 t 
мії Ө, = Ө апа Ө, = Ө + 180°. Substitute Roark’s o(@) into 


the above equation gives: 


14 . 
a + =sin0 
EE AS w'rCa -"sinO 5 t 

г а + sino, 

4 
‚ а + zsinO. 
+ Pr (2 m rca + sino, >t 
г а + sinó, 


Cancelling some terms shows the equality: 


да = 2a Q.E.D. 
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This must have been how Roark arrived at these relationships. 


This is how the solution will be checked. 


THE SOLUTION 


The method outlined above requires the integration of 
several quantities which are not readily found in a table of 
integrals. The integration was done using Simpsons Rule 
(231. It was found, using Known integrals, that 50 steps 
with in the interval O to 2n gave good accuracy (5 to 6 
significant digits accuracy). A program was written to 
accomplish the integration and the subsequent matrix 
manipulations to solve for {u} and tw). 


- 


To check the output from the program, a spreadsheet was 








set up to handle: 
WO = 5 м $ ont O? pega Ww $ 
К 50 
== ulo) _ - E 
uo -2u $6 P zu $ 


The spread sheet also calculated Се Еф» O9» and Оф: Аз а 
check, 0% was averaged and used to calculate a pressure to 
compare with the input pressure. 
THE RESULTS 

The results are disheartening. In the sample output 


provided, only one of the stress values used as a comparison 





came close to the expected value. 0% should have been 
constant and have a value of 5000 psi. It was not constant 
but it did manage to equlibrate most of the applied pressure 
€99.9%). This is pretty good for an approximation. The JG 
is another story. 


da is too low and does not follow the solution in Roark: 


4: 
DET Р а + ;sin® 
t a * sinO 


Фа 


The following are the default values in the program: 


GEOMETRY 

Radius of rotation : 8 inches 

Radius oft the circle : 2 inches 

Thickness of the shell : .02 inches 

MATERIAL 

Poisson’s ratio : .3 Modulous of elasticity 
3E+07 

SOLUTION 


Number of steps in the Simpsons 


integral : 50 Number of terms to be 
used in 

the series approximation : 10 
Pressure :-100 psi 
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The program output for м. апа u_ are 

m ы” u 

1 -1.213117Е-03 -1.137496E-03 
2 6.043412E-09 1.491326E-09 
3 -1. 034239E-04 -1.151281E-05 
4 3.250964 E-09 2.026468E-10 
5 =3. 933767E-05 -1.579169E-06 
6 3. 328626E-09 9.22/051E-11 
% -1216143E-05 -2.382749Е-07 
8 3.831875Е-09 5. 9864E-11 

9 -2.034983Е-06 -2.535156Е-08 
10 2.81034Е-09 2.817199Е-11 


Notice that the even numbered terms, which are to be set 


to zero as indicated earlier, are several orderes of 


magnitude below most of the odd values. 


The reconstructed output from the program is as 


follows: 





О 
15 


ЗО = 


45 
60 
(44 
90 
105 
120 
135 


165 


180 - 


195 
210 
225 
240 
255 
270 
285 
300 
315 
ззо 
345 
360 


Мое: 


The spreadsheet out put to check the solution is: 


~) A 


ho MM O he BP hehehe 0 


EN 


1? 
Е 
1-ъ 





с OON 


y 


Cin» 


. OE*00 
74:04 
sd 
. OE-04 
10Е-03 
.1Е-03 
заЦесаше: 
ШЕ ОС 
-0E-03 
. 0Е-04 
150 -7. 


2Е-04 


.4E-04 
. ÓE-08 
.4E-04 
“2Е-04 
. OE-04 
. OE-03 
=1E>03 
-TE-03 
1Е-03 
. OE-03 
. DE-04 
.2Е-04 
.4Е-04 
52Е-07 


1 
POOR PROQOO IE NN 


ою 


-100.00 
0.30 


3.3E+07 
8.00 


2.00 
0.02 


Ow 

50 

in 
Eu 


. BE-03 
.4E-03 
.1E-04 
. TE-04 
. 2E=04 
.2Е-04 
-ЗЕ=06 
. ZE-04 
. ZE-04 
. f E-04 
.1E-04 
.4E-03 
. SE-03 
.4 E-03 
.1E-04 
. VE-04 
.2E-04 
.2Е-04 
. BE-08 
. 2E-04 
.2E-04 
>СЕ- 04 
.1Е-04 
.4Е-03 
. BE-03 


psi 
E 
in a 
in 
in 
u 
Cin? 
-1.2Е-03 


.1E-03 
.8Е-04 
. BE-04 
.4E-04 
. 8E-04 
. BE-08 
. BE-04 
.4E-04 
. BE-04 
. BE-04 
.1E-03 
. 2E-03 
.1E-03 
. BE-04 
. BE-04 
. 4E-04 
. BE-04 
.4E-08 
. BE-04 
.4E-04 
. BE-04 
. BE-04 
.1E-03 
.2Е-03 


Due to the precision of 


O ROO CER BRB O DA PRO 


-99.90 < Sln Press Cpsid 


3E+07 


4. 


Su 

ФӨ 

іп 
zc 


. OE*00 
. 2Е-04 
.8E-04 
.4E-04 
. 6E-04 
. OE-03 
.1E-03 
. OE-03 
. 6E-04 
.4E-04 
.8Е-04 
.2Е-04 
.2Е-08 
.2E-04 
.8Е-04 
.4Е-04 
.6Е-04 
. OE-03 
.1E-03 
. OE-03 
. 6E-04 
.4E-04 
. BE-04 
. 2E-04 
.8Е-07 


Фе 
Cpsid 


-1461 
-1710 
-1946 
-2123 
-2258 
-2345 
-2375 
-2345 
-2258 
-2123 
-1946 
-1710 
-1461 
-1269 
-1085 
-894 
-747 
-659 
-630 
-659 
- 747 
-894 
-1085 
-1269 
-1461 


99 
Cpsid 


-4870 
-4738 
-4607 
-4405 
-4248 
-4159 
-4131 
-4159 
-4248 
-4405 
-4607 
-4738 
-4870 
-5190 
-5498 
-5654 
-5769 
-5854 
-5886 
-5854 
-5769 
-5654 
-5498 
-5190 
-4870 


the machine being used, the 





forced compatability conditions appear to be non-zero. The 
forced conditions at © = 900 and 270° are significantly less 


than the rest of the solution and therefore satisfactory. 
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CHAPTER 6 
CONCLUSIONS 
The apparent simplicity of the toroidal shell belies the 
complexity of the mathematics to analyse this shell 
structure. The analysis of simpler shells Cplates, 
cylinders, and spheres) has the advantage of the geometric 
curvature remaining either constant or at lease maintaining 
its sign positive or negative. The geometric curvature of 
the toroid changes from positive to negative as © goes from O 
to 21. 
The impact in the change of sign in the radius 
of curvature about the Z axis has been seen in the results 
of this work and in applying other*s solutions. This impact 
is emperical and could be the subject for further 
investigation. The first blatant statement was found in Tsui 
{24} where he states, "...the relevent differential equations 
to assume singular solutions, and consequently, problems 
involving the crowns such as a complete toroid cannot be 


solved." The culpret was writing the differential equation 
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in terms of «XTsui's nomenclature»: 


15 b C1 + A 2 
sino 





(This papers nomenclature): 


R 


ре г + C=) “sin®) 


1 
sino 


To avoid the term, Tsui shifts variables - both 





dependant and independant. - and provides influence 
coefficient matrices. (These are outputs from numerical 
analysis. > 

As mentioned earlier, Tsui addresses the significant 
differences between the partial toroidal shape and the 
complete toroid. Others may have addressed partial toroidal 
shapes tangentially but’ Tsui once again is very blunt in his 
statements. Again remember Tsui is intending his work as a 
practical handbook. His intended audiance wants practical 
answers to real problems. Others, Flugge [25] 
specifically, attempt to give solutions but he adds, "...this 
solution cannot be realized...because it again leads to an 
incompatability of deformations..." 
CURVATURE 


In attempting to point the finger at why the solution 


is so difficult, the answer comes back to curvature. This 
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work has seperated curvature into geometric (unloaded> and 
loaded curvature. This writer feels that the geometric 
curvature is the culpret. Notice the expression for the 


radius of curvature in the Фф direction given in chapter one: 


R 
sino 


For the range 0% < e < 180%; sinó > O and p > 0. Conversely, 





+ г 


for 1809 « e « 360°; sine < 0 апа р < 0. At © = 0° and 180° 

sino = 0 and pç = + о. Two things were observed. 
A ln previous solutions obtained in this study- ie. 
worse than the one presented - the stress outputs for 
the two regions differed markedly between 0% <e< 
180° and 180% < e < 360%. This could have 
been due to a lot of reasons, but this was also 
coincident with 
В. Flugge’s solution for the u displacement [26] 
was being investigated. This tended to blow up from 
Ө = 180? to 3609709 due to a term ln(tan$>. Other 
problem terms existed Ccot®) which agrivated the 
solution at the points © = 0° and 180°. In just 
“playing with the solution, the points of curvature 


change kept comming back as critical points requiring 





further attention. 

Something happened at O = 09 and 1809. In a cylinder 
these two points are no trouble - but the radii of curvature 
are « and r - and are constant! In a sphere these two points 
are no problem - but the radius of curvature is r in all 
directions - and constant! In the toroid the radii of 
curvature at both points are p = + w and r. The key 
difference is that p is changing from p > O to p < O and at 


ф Cor > just happens to be passing through 


© = 0° and 180°, x 
zero. 

It was mentioned earlier in Chapter 4 that it required 
two reference points to define the toroid. Some individuals 
may take exception to that statement, but here is where it 
comes into play. Another way to look at p, the radius of 
curvature in the $ direction, is that p is the length of the 
arm joining the point in question Csay point Pd and the axis 
of symmetry (the Z axis in our case) and coincident with the 
local unit normal vector at point P. The angle © can be 
measured either at the axis of rotation of the small circle 


forming the shell, or at the intersection of the line p and 


the axis of symmetry for the entire toroid. Now when © = 50° 





for instance, the geometry looks like that shown in Figure 
16. As © decreases toward 0, the end of p must travel 
farther and farther down the axis of symmetry. All this time 
p is getting larger and xp is getting smaller. When © gets 
to zero plus (0°) the situation is one of approximations. p 
is the hypotenuse of the right triangle with the right-angle 
at the origin. Sin® = х from simple trigonometry, but аз р 
goes to a, sin® goes to zero and hence © must go to zero 
also. 

Continuing © around to 0 , the closest point on the axis 
of symmetry in line with the normal vector is now at +œ. 
This now gives p a direction opposite to that of the local 
normal vector, hence a minus sign, and as © goes further 
negative, p is how smaller in magnitude. CSimilar arguments 
can be made for the bottom crown where © = 180% The sole 
reason for the + w trip in o, and hence the sign change in 
x is the offset R. Somehow this must be related to the 
difficulty in solving the complete toroid. 
GEOMETRY 

The other subtility of the toroid is the way the 


mathematical representation of the surface impacted on the 
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attempt to integrate for energy over the toroid surface. 
This work specifically avoided some problems encountered by 
other investigators, significantly the expression for the 
curvature in the o direction. The early development in this 


work specifically avoided the term С >, selecting instead 





Sino 


an expression containing the term Са + sin®. Force balance 
and constitutive equations were developed using the latter 
expression. When the constitutive equations were plugged 
into the energy expression CV), no singularities existed. 
This was good news. 

Unfortunately, the Ca + sin®) term also appeared in the 
expression for the differential surface area CdS). The bad 
news was that. many of the elements in V contained 
combinations of trigonometric functions not readily found in 
the CRC Math Tables [27]. This is what forced the numerical 
integration of the elements in V. 

This heavy dependence of the trigonometric functions in 
describing the problem has an astounding impact on the 
selection of an assumed functional for the energy method 
solution. Any combination of CA 8 sin® ог ce” sine goes to 


zero when integrating from O to 2n. Remember that the dS 





term contains a sinO term so all integrals contain at least 
one Sin® term. 

It should be pointed out that the solution to a 
symmetric problem can be obtained by solving only one of the 
unique parts of the problem. This was not done in this case 
because it was hoped that the program being developed would 
be more general in nature, actually leading to orthotropic 
analysis of a stiffened toroid. 

The only redeeming guality of the integrals also comes 
from geometry. Geometry terms containing cos®, cos*®, sin@, 
sin*@, and «а + <1п) and thier products appear in numerators 
and denominators through out the integrals. Mister Simpson 
and the digital computer were allowed to seperate the wheat 
from the chaff. 

An attempt was made to develope a non-trigonometric 


functional which would satisfy the conditions of 


connectivity stated earlier and to meet connectivity at 
© - 0973609. The conditions to be met were: 
fCO> = FC360) 


г (0) + Г (360) 


ГС902 0 
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f (900 =0 
fC270 = 0 
f (2702 = 0 
This resulted in eight conditions. A functional of the 
following form was generated: 
f(m6> = A + B(m6) + CCm6>* + DCme>* 
+ ECm6>* + FCm6>* + GCmO>* + H<me>" 
A was arbitrarily set to 1 (to force a non zero solution at 
o 0°). Unfortunately, by the time fCmq> was processed, 
the functional looked the same for all values of m. The 
coefficients all changed as m changed, but the shape of the 
functional remained unchanged from one value of m to the 
next because of the eight conditions imposed upon the 
functional. This lost the variety of the solution that one 
expects to obtain using a varying m. See Figure 17. 
It is obvious that more complex functionals are needed. 
This work included an initial evaluation of several 
funct.ions for this task. A considerable amount of time was 
put into investigating the hypergeometric functions, 
particularly the Legendre functions. In addition to serving 


as a functional for u and w, it was also hoped that this 


— 0 = 





function would satisfy the form of the differential equation 
for the toroid and hence yield an exact solution. (It was 
partially the anticipation of using such a complex function 
that drove the layout of the program (see Appendix) to have a 
seperate generating section for the functional.» The 
Legendre function, or any other hypergeometric function, was 
not used due to compatability conditions being harder to 
achieve. The next step in this progression would be the 
assumption of a combined trigonometric function for u and w, 


Cie. An sincm@) + B cosCm@) >. 
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Figure 16 


Radius of Curvature 


Ди, 








Figure 17 


Generated Functional 
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CLOSING 
In Chapter 1 it was stated that the geometric curvature 


and the impact of shell geometry on the assumed functional 
were the two problem areas in this thesis. Of the two, the 
former is the more significant in this writers opinion. The 
analysis behind Figure 15 is crude, but it is an attempt to 
represent in graphics and mathematics what the shell is 
doing. Being unable to do this is what makes the solution so 
difficult. 

In a dynamics class, one student asked the professor if 
another mathematical technique could be used to analyse a 
spring. The professor thought briefly, then answered, "You 
can do all the math you want but you have to think like a 
spring. The spring Knows what it is doing.” This writer 
feels a little short on thinking like a toroid, but also 


feels in good company. 


a= 





APPENDIX 


THE PROGRAM 


This program is to set up the geometry, the energy. 


expression, Simpson’s multipliers. and perform the matrix 


manipulations to solve for 


he displacements £u? and tw». 


20 due UE PRID: INPUT "WHAT OUTPUT FILE NAME ';A8 
“ 
1000 NN= 50 *This is the number of steps to be used in the 
Simpson's integral | 
1010 ММ=10 *This is the number of terms in the series for 
displacements u ans wv. 
1020 NU=.S *Poisson’s ratio 
1030 R=2 XR from the text. 
1040 A-4 *a from the text. 
1050 H=.02 жі from the text. 
60 E=3E+0/  *Young's modulus 
P=-100 pressure in psi. 
GESTERINT: PRINT: ERINT 


DOSE casa 


ON UT IE CO NOH ОООО-ЛОХҚООМНОООЗООООМНСОО OOM 


БА БА Ба араа А Қ A a mE А |» 
me Gee Оо СС) сс ОЮМ ммм ѓе неее OOOO 


CC) OC -есоссесссссссссоа 


СЛО О ЕСО 0351 


PRINT “The following are the default values in the 


rogram: 
PRINT: PRINT 
PRINT “GEOMETRY” 


PRINT “Radius of rotation >: AKR:" inches” 


PRINT "Radius of the circle : ;R; inches" 
PRINT "Thickness of the shell : ;H;" inches” 
PRINT: PRINT 

PERINI “MATERIAL” 

PRINT "Poisson's ratio 0: NU 

PRINT "Modulous of elasticity WE 

PRINT: PRINT 

PRINT "SOLUTION" | 

PRINT "Number of steps in the Simpsons" 

FEINT ” integral 2's NN 

PRINT “Number of terms to be used in” 

PRINT " the series approximation MH | 
PRINT "Pressure "sP;" psi" 


FRENI: PRINT 

INPUT "Do you wish to change any of these (Y/N): ': 
Е аео НЕХ 1520 

PRINT: PRINT 

РЕТИ “GEOMETRY” | 

INPUT "Radius of rotation Cinches) :“, ВК 

INPUT "Radius of the circle (inches) : "R 


A=RR/R 
THEN 1380 


S$ 


Im] 
N OF ROTATION MUST BE GREATER THAN THE 
A с" 


RADIUS 
18002077 THE CTRCLE. :GOTO 1300 


INPUT "Thickness of the shell Cinches) :";H 
PRINT: PRINT 

РЕШИ "NATEKIAL" 

INPUT "Poisson s ratio a NU 
INPUT “Modulus of elasticit AE 


PRINT:PRINT1440 PRINT "SOLUTION" 
PRINT "Number of steps in the Simpsons" 


m Aue 





нн E ы 
CHINO) N UI eee 
мМ -O09-10 
oO “See eae 


~Ja de go 
COOOC 


the partial of V with respect to u, 


representing the energy due to S бак term. 
scratch matrix. 
matrix manipulation subroutines. 


INRUT ~ integral S NN 
PRINT "Number of terms to be used in" 

INPUT the series approximation MM 
INPUT "Pressure (psi) Zr 


"THIS SECTION IS TO SET UP THE PARAMETERS FOR THE 
„INTEGRATION SECTION 


DT=6. 283185307 #/NN 
I=NN+1 

DIR PCI. 3, MM) 

DIM JCI) 
Arrays MD, ME, MF, 


*Defines the functions ф. 
*Simpson?s multiplier 


and MG ко to the matrixes in 
and w MP is the matrix 


MS is a 
is a series of matrixes used in the 
MW and MU are the answers. 


Mxr 


1580 DIM MDCMM,MM>: DIM MECMM, D DIM МЕСММ, ММ» 

1590 DIM MGCMM.MM):DIM MSCMM.MM) 

1600 DIM MX=CMM, MM>: DIM MYŻCHM, ММ»: DIM MZ4CMM, MM) 

: DIM МІ#СММ, MM) 

1610 DIM MPCMM): DIM MWCMMD: DIN MUCMM) 

1620 CLS: PRINT: PRINT: PRINT "GENERATING FUNCTIONS" 

1630 FOR I = 1 TO NN+1 

1640 T=CI-1)*DT 

1650 FOR M = 1 TO MM | 

1660 ЕСІГІ M»zsINCMM "This is the function 

1670 PCI.2,MD=M&COSCM&TD *This is the first 
derivative 

1680 PCI.3,M0=-M 2%5INCM*T) "This is the second 
derivative 

1690 NT A 

1700 ДЕ 

1710 Б PE INTCI/Z2»5»0 THEN JCI»22 

S. 0 NEXT I 

1730 теке 1: JONN+1)=1 


Integrate using Simpson’s Rule 


2000 
2010 
2020 


FOR I 
0 


> INTEGRATE FROM O TO 2*PI 
ІС-3.14159ж%СЕжН/(1-У07222жрТ/У3 
BC=H* 2/C6*R* 3) 
PC=-2*3.14159%*P*R* 2*DT/3 

= 1 TO NN+1 

T=C 1-1 3xDT 

S=SINCT) 

C=COSCT) 

CLS: PRINT: PRINT: PRINT INTEGRATING ON ЭЛЕ 551 
ЕИ ANGLE *.Ta57.296,  DEGREES 
AA=A+S2102 HH=H*2/C12*%R*2) 
CC=C/AA 

CA=2*NU*CC 

CB=CC” 2 


CD=2*NU*S/AA 
AX=2*S*C/AA’ 2? 





2130 FOR K = 1 TO MM 
2140 MPCK=MPCK)-JC1>*PCI,1,K)*AA 
2150 SQ=1/AA 2 
2160 TR=2*NU/AA 
2170 BX=H*3*C1+2*NU*C/AA+(C/AA)*2)/C12"R*2) 
2180 FORM = 1 TO MM 
2190  MDCK,M)=MDCK, ОША A ECL IR 3, O*PCI,3, MD 
AXPCI.3,KO*PCI,2,MÓ ' 
*2xPCI 2 ,KO%PCI 2.MD)J*AA 
2200 MECK,M)=MECK, + CISC KPC T 1,MxPc1 3 
+САЖСРСТ.З.ЕУжРС1.1, МУже 
+РСТ ЫН 2 KXC) 
АБ > opel 1 
PLA Морт 
*CAX( Ра 3 Ы 5. M) 
+РСТ 3 МУжРСТ 2,K00 
| | +2жСВЖЕСТ 2, Ова. 2,M)JJ*AA 
2210 MFCK. M)=MPCK, M+ СТ) ЖС2ЖРСТ,З3. МУжРСТ.1.КУ 
£XXNUXCPCI. Трето 3. 552 
+PCI 2 Mee CH 1 ЮЖОУАА 
+АХЖРСТ. 2. МУЖРСТ.1 
к ыы ATTE) 
Canc PC 3 AC 2 
£PCI.3, KÓWPCI 
2xGBHPCI, 2 Be. 2. K))>*AA 
2220  MGCK,M)=MGCK,M)+JCI)*C2*PCI я *PCI IR 
+2xCDxPCI 1 MOKPCT >” 
жест. 1. T3268 1 ок: L4 AA" 2 
THHX*C2XPCI,3.MO*PCI.3 
+CA*(PCI З Web, 2, KD 
+PCI.3.K)*PCI 2 
анын ы 2. ҚОУОЖАА 
2230 МЕХТ М 
2240 NEXT K 
2250 NEXT I 
2260 ? 
3250 > 
2280 FOR K = 1 TO MM 
2290 MPCK)=PC*MPCK) 
2300 FOR M = 1 TO MM 
2310 MDCK,M)=IC*MDCK, MD 
2320 MECK,M)=IC*MECK, MD 
2330 MFCK.M)=IC*MFCK.M) 
2340 MGCK.M)=IC*MGCk.M) 
2350 NEXT M 
2360 NEXT K 
3000 *SOLVE THE MATRICES FOR w_ AND u_ 
3010 ( CLS: PRINT: PRINT: PRINT "SOLVING" 
3030, FOR k = 1 TO MM:FOR M = 1 TO MM: MX4CK, M)=MDCK,M): NEXT 
3040 GOSUB 10000 " жжжжж Invert MD  ebiok 
3050 FOR K = 1 TO MM:FOR M = 1 TO MM 
3060 MY#CK,M)=MECK,.M):MX#CK, MD=-MIACK, MD 
3070 NEXT M: NEXT E 
3080 GOSUB 11000 ° жжюжж = C-MD 7-1) #ME soloiok 
3090 ҒОР К = 1 TO MM:FOR M = 1 TO MM 
3100 MY#ZCK, M>=MZ# CK, M): MX#CK,M)=MFCK, M) 
3110 MSCK, MD=MZ#CK, Ñ) жююк Saving this for Um 
solution кок 
3120 NEXT M: NEXT K 
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GOSUB 11000 ' жк МЕЖС-МО7-1ОЖМЕ olelok 

FOR K = 1 TO MM:FOR M = 1 TO MM 

MY#CK, MD=MGCK, MD: MX#CK, MD=MZ#CK, MD 

NEXT M: NEXT K 

GOSUB 12000 ? ЖЖЖ MG — MFxC-MD"-10%ME — Wetotot 
FOR K = 1 TO MM:FOR М = 1 TO MM: MX#CK, M)=MZ#CK, MD 

_ :NEXT M: NEXT K 


"ROV REDUCTION 
iride is MX*CK,M> and MPCK) 
2 ӨЛІ ut is the altered MX#CK,M) and МРСКО 
FOR 1 TO MM-1 
FOR K = S+1 TO MM 
=-MX#CK.D/MXZCS. DD 
PCK)= =O*MPCS)+MPCK) 


FOR M = 1 TO MM 
MX#CK. MO=Q*UMX¥CS, MO+MX¥CK, MD 
NEXT M 
NEXT К 
АБАТ S 


SOLVE FOR MWCM) 
Continuing dw using MA*CKk,M> and МРСК? from before 
FOR S = MM TO STEP -1 

MWCSD=MPCS)D 

FOR ^ = S+i TO MM 

MWCSD=MWCS)—-MWOMD *#MAX FCS. MD 
NEXT M 
MWCSD=MWCSD/MAFCS, SD 
NEXT 


2 
= 
> 


ЖЕСЕУЕ КОЕ МОСК? 

FOR k = 1 TO MM 

КОК М = 1 ТО ММ 
MUCK>=MUCK>+M=CK . MO*MVCM) 
NEAT И 


NEXT k 
i PRINT OUT RESULTS 
PRINT 7," Mm ",* Um " 
FOR M = 1 TO MM 
BEDNT. M. MWCM5. MICMD: PRINT 
МЕХТ М 
ОРЕМ "0", #1, АЗ E section puts the output onto 
isk. 
PRINT £1,"NN 5» "МА." MM = ';MM 
PRINT #1. "RADIUS ="; AR, radius =";R 
PRINT +1, 
PRINT #1. ova Um ” 
FOR M = 1 o MM 
КЕЧЕЕ М ПАСМО MUCHO: PRINT 
NEXT M 
CLOSE 
END 





This section contains routines used in the main program. 


060 


© 
O 


UE WHF OOD ONOQUF WHR O 


fa ja js js js jas fas p as as fa fas a jas jas jas as a pajas jas jas jas jas as pas ps 
SIGN 
е» «>» зл] а» | <> «жь «=» > е «в а» | [ж ез | «т | «к» | нэ 


ajajaja aja ao ao Ссс О с == оссо 
DOBRO DO RID DODO ROR RRP PRP eRe Pr OO 


М) 
CO 
O 


O 
vr 
сс 


МММ PRERERRERRRRRRR ООО 
oo тасос сее 
HO MOON NO MC 
сто =) 6 ОСЕ СС 


N 

об 
WN 
OO 


NNN 
OOO 
ON Ub e 
ооо 


jà pà pd pd pd pd pd pd pd Һә Һә Һә papas 
N 
O 
-] 
O 


“input iS МХЕ 
200 utput i is MICK. M) 
TO MM 


"Mak Ke 1 
NEXT k 


T 
FOR S = 1 TO MM-1 
FOR K = S+1 TO MM 
О--МХФСК, ¿SIÓMXECS. 
FOR M = i TO MM 
MX#CK, м» =O*MX#CS, 
ЕЕЕ MJ=Q*MIECS, 


" INVERSION m A ие SS 


(и) 


) 


= 
VY 


+MXFCK, M) 
«МІРСК,МО 


MM TO 1 STEP -1 

= S-1 TO 1 STEP -1 

Us Ny K. E)/MX#CS S) 
OR M=MM TO 1 STEP -1 
MN SCK. MD EQXMXECS, D *MXECK, ND 
MIZSCK, D zOXMISCS. ND&MIZCK. MD 
NEXT M 

NEXT K 

NEXT S 

FOR K = 1 TO MM 

FOR M = 1 TO MM 
MI#CK,M=MI#CK.M/NX#CK,K) 

NEXT M 

NEXT K 

RETURN 


"MULTIPLY TWO MATRICES 

"Multiply MX#CK.M) times MY¥CK,M) in that order! 
"The output is MZeck,M) 
FOR K = 1 TO MM 

FOR M = 1 a MM11050 MZ+CK,M>=0 

FOR S TO MM 

MZ=CK. "n MZECK MO+MXFCK SIO*MYŚCS MV 

ДЕ ABSCHZECK, M3><.0000005 THEN MZ#CK,M)=0 


RETURN 


"ADD TVO MATRICIES 

"Add MX#CK.M) to MY#CK,M) 
"Output is MZĘCK, M) 

FOR K = 1 TO MM 

FOR M = 1 TO MM 

MZeCkK.MO=MX/CK, MO+MYŚCK, M) 

NEXT M 

NEXT K 

RETURN 
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